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ABSTRACT 

High-quality integrated radio profiles of some pulsars contain bifurcated, highly sym- 
metric emission components (BECs). They are observed when our line of sight tra- 
verses through a split-fan shaped emission beam. It is shown that for oblique cuts 
through such a beam, the features appear asymmetric at nearly all frequencies, except 
from a single 'frequency of symmetry' Vsyim at which both peaks in the BEC have the 
same height. Around Vsym the ratio of flux in the two peaks of a BEC evolves in a 
way resembling the multifrequency behaviour of J1012+5307. Because of the inherent 
asymmetry resulting from the oblique traverse of sightline, each minimum in double 
notches can be modelled independently. Such a composed model reproduces the double 
notches of B1929-f 10 if the fitted function is the microscopic beam of curvature radi- 
ation in the orthogonal polarisation mode. These results confirm our view that some 
of the double components in radio pulsar profiles directly reveal the microscopic na- 
ture of the emitted radiation beam as the microbeam of curvature radiation polarised 
orthogonally to the trajectory of electrons. 

Key words: pulsars: general - pulsars: individual: J1012-I-5307 - B1929-I-10 - Radi- 
ation mechanisms: non-thermal. 



1 INTRODUCTION 

Double absorption features have so far been observed in in- 
tegrated radio profiles of B 1929+ 10 (Rankin & Rathnas- 
ree 1997), J0437-4715 (Navarro et al. 1997) and B0950+08 
(McLaughlin & Rankin 2004). They are the peculiar 'W- 
shaped features observed in highly polarised emission. The 
minima of the 'W approach each other with increasing fre- 
quency v. They have a large depth of 20 — 50% (Perry & 
Lyne 1985; Rankin & Rathnasree 1997). 

Initial attempts tried to explain the features in terms 
of a double eclipse, with the doubleness generated by spe- 
cial relativistic effects (Wright 2004; Dyks, Frackowiak, 
Slowikowska, et al. 2005). Dyks, Rudak & Rankin (2007, 
hereafter DRR07) proposed that the feature represents the 
shape of a microscopic beam of emitted coherent radiation. 
It has been suggested that this beam is also observable in 
pulsar profiles as a double emission feature, hereafter called 
a bifurcated emission component (BEC). Examples of such 
components can be found in the pulse profiles of J0437— 4715 
and J1012-I-5307 (Dyks, Rudak & Demorest 2010, hereafter 
DRDIO). However, the beam that had been proposed ini- 
tially by DRR07 - the hollow cone of emission caused by par- 
allel acceleration - could not reproduce the observed large 
depth of double notches (20 - 50%). 

DRDIO have shown that the drop of flux in the minima 
of double notches by several tens percent imposes constraints 



on the general structure of the emission beam. The flux drop 
in the 'W must be caused by either a small non-emitting 
region (wherein the coherent amplification fails), or by some 
obscuring/intervening/eclipsing object. For the hollow cone 
beam, the decrease of fiux expected in such a model is max- 
imally a few percent (see figs. 2c, and 4 in DRR07). This is 
because within the hollow-cone beam a considerable amount 
of radiation is emitted along many non-parallel directions. 
A single line of sight is then capable of receiving consider- 
able contributions from sideway emissions. For this reason 
a laterally extended, two- (or three-) dimensional emission 
region contains many places that backlit the notches and 
make them shallow. Numerical simulations have shown that 
this "depth problem" cannot be naturally overcome through 
manipulations with geometry of the non-emitting (or ab- 
sorbing) region (see fig. 5 in DRR07 for sample results). 

In DRDIO it has been recognized that the depth ap- 
proaching ~ 50% results from the fact that the emission 
is intrinsically two-directional, or, that the beam is double- 
lobed. In such a case, at any pulse longitude the observer 
simultaneously detects radiation from two separate places 
located somewhere within the extended emission region. At 
the phase of minimum in the notches, only one of the two 
places is eclipsed (or is non-emitting), hence the 50% upper 
limit for the depth of notches. 

When the double-lobed beam is carried along with the 
electron, a split-fan beam is created, as shown in Fig. [1] 
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Figure 1. The principle of creating double features in radio pulse 
profiles: the line of sight traverses through a split-fan emission 
pattern, produced by outflowing particles that emit a double- 
lobed beam. For simplicity, the figure shows only a single trajec- 
tory of electrons flowing in a stream of flnitc width. 

In this paper we show that when the line of sight cuts 
through the stream at an oblique angle, a specific change of 
the EEC's shape is expected across the frequency spectrum 
(Section 2) . The resulting asymmetry of the BEC allows for 
a piece-wise fitting of double notches with a physical beam- 
shape model (Section 3). 

One important difference between the stream-cut origin 
of components and the conal view of profile formation is the 
following: the altitude of detectable emission corresponds to 
that part of the stream, which is nearly tangential to the line 
of sight, whereas in the conal model, the altitude is fixed by 
a preselected value of v (the tangency condition determines 
the azimuth of detectable region, but not the altitude). In 
the case of the stream-cut scenario, when the line of sight is 
not in the plane of the radio-emitting stream, no (or little) 
radiation is detectable, irrespective of v. When the sightline 
crosses the stream plane, only a localized piece of the stream 
is detectable, namely, the part that happens to be (nearly) 
tangent to the line of sight. Thus, the 'observable' altitude is 
frequency-independent, because the condition for tangency 
does not depend on the frequency of observation. Conse- 
quently, if the components observed in pulse profiles arise 
from stream cutting (whether they appear double, or sin- 
gle because of convolution effects), they will not exhibit the 
radius-to-frequency mapping (RFM). That is, the compo- 
nents will remain at fixed, frequency-independent locations 
in pulse profiles observed at different v. 

Such a stream-cut origin of components allows us then 
to explain the lack of RFM in millisecond pulsars (MSPs; 
Kramer et al. 1999). Stream-cut viewing geometry is more 
likely in MSPs, because the polar cap of MSPs is much larger 
than for normal pulsars. Therefore, any transverse motion of 
limited magnitude, like the ExB drift, is less likely to create 
the hollow-cone emission pattern. Moreover, magnetic field 
lines in the MSPs flare strongly away from the dipole axis, 
which facilitates cutting of sightline through streams. 



2 OBLIQUE CUT THROUGH SPLIT-FAN 
BEAMS 

The split-fan emission of a stream, when projected on the 
sky, creates an elongated double track, such as the slanted 
pattern shown in Fig. (2] There, the sky-projected emission 
is presented at three frequencies: vi (two left contours: LI, 



and Tl), 1^2 (grey central contours: L2, T2), and 1/3 (two 
right-hand side contours: L3, T3). For clarity of the figure, 
the contours mark only the brightest (topmost) parts of the 
emission pattern. There is actually considerable (and de- 
tectable) level of emission extending all the way across the 
figure at each frequency. Passage of sightline through the 
pattern, eg. along the horizontal line in Fig. (2] results in es- 
sentially double shape of observed component, regardless of 
the frequency of observation (see the triple inset). Assum- 
ing some residual RFM, for decreasing f the double track 
has a maximum at increasing distances from the polar cap, 
ie. the location of the peak fiux moves from right to left 
in Fig. [2I For appropriate traverse of sightline, there ex- 
ists a frequency (denoted 1^2 in Fig. [21 and hereafter called 
the frequency of symmetry, i^sym), at which both maxima 
in the BEC have the same peak flux. Note that for non- 
orthogonal cuts (5cut 7^ 90°) through the split pattern, the 
line of sight does not cross through the maxima of the t'sym" 
pattern (marked 'L2' and 'T2' in the figure). Even at this 
frequency of symmetry the sightline misses the leading max- 
imum L2, moving on the 'far side' of it. That is, the sightline 
passes through the stream at a larger radial distance ta than 
that of the peak emission: r™'"'(i/2) = r^'^^ < va- The trail- 
ing maximum of the i/g-pattern is passed- by at vb < r™'"', 
ie. below the radial distance of peak flux at this frequency. 
Thus, even when peaks of an observed BEC have perfectly the 
same-height, they nevertheless correspond to different emis- 
sion altitudes (for Scut 7^ 0). For this reason, despite having 
equally high peaks at Vsym, the observed BEC will not be 
ideally symmetric: the peaks can have somewhat different 
widths due to a difference in curvature radii (see below). 

At a lower frequency 1^1 < Vsym the BEC will have the 
leading peak brighter than the trailing one, because our line 
of sight passes closer to the leading maximum in the 
pattern (LI). Using the markings of Fig. [51 at frequency 
vi < Vsyui we have |L1 A| < |T1 B|. This holds (ie. the lead- 
ing peak is the brighter one at 'subsymmetric' frequency) if 
the polar cap (region of lower altitudes) is located on the 
right-hand side (trailing side) of the observed feature, as in 
Fig. [2I That is, the leading peak is expected to dominate 
at v < i^sym if the BEC is observed as a precursor; for a 
'postcursor' position of the BEC, the opposite is expected: 
a brighter trailing peak a.t v < I'sym. 

Conversely, the trailing peak of a BEC becomes brighter 
at V3 > Usym, because |A L3| > |B T3|, see Fig. [H Thus, 
when the BEC is inspected at increasing frequencies, the 
dominating (brighter) peak switches from the leading to trail- 
ing position at Vsym (see the inset of Fig. [2} . Observations of 
Kramer et al. 1999 (fig. 5 therein) provide evidence of such 
behaviour in J1012-I-5307. For the precursor location of the 
BEC, it is the leading peak which is expected to be brighter at 
low V, again in consistency with the profile of J1012-I-5307. 

2.1 Asymmetry of the width at Vsym 

At the "frequency of symmetry" the maxima of BEC have 
precisely equal height, but the symmetry is approximate, 
because their widths are slightly different. The reason for 
this is that the width-scale of the BEC also depends on the 
curvature of trajectory of electrons flowing in the stream (in 
addition to the dependence on 5cut). Low in the dipolar held, 
the curvature radius p increases with altitude, as a result 
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Figure 2. Origin of asymmetry inversion in bifurcated components of pulsar profiles. The three pairs of contours present sky-patterns 
of flux at three frequencies ui < Usym < ^3- The patterns result from emission of a stream extending quasi-diagonally towards the 
bottom-left corner of the graph. Observer's line of sight crosses the stream along the horizontal line. For more details see Section [2] 



of which the emitted beam becomes narrower. The leading 
half of the BEC originates from larger altitudes (see Fig. ^ 
and should therefore be narrower than the trailing part, for 
which the altitude and p are smaller. Assuming the oblique 
cut and the normal RFM, the bifurcated precursors are then 
expected to have the leading wing narrower than the trailing 
one. Weak traces of such asymmetry can be discerned in the 
BEC of J1012+5307 (see fig. 4 in DRDIO). 

This inference has been verified with a numerical code 
that follows beam-resolved emission of electrons moving 
along narrow streams in rotating dipolar magnetosphere. In 
the numerical calculations illustrated in Figs. 3 and 4, we 
assume that the B-field has the geometry of the rotating 
retarded vacuum dipole as given by equations Al - A3 in 
Dyks & Harding (20040- Thus, we do not take into account 
the modifications of B, that are expected for plasma-filled 
magnetosphere (see eg. Bai & Spitkovsky 2010 for the force- 
free case) . The usual static-shape dipole is assumed only for 

^ We take the opportunity here to warn that fig. 1 in Dyks 
& Harding (2004), that compares the retarded and static-shape 
dipole, was corrupted by errors in plotting routine and should be 
dismissed. All formulae in that paper, however, are perfectly free 
from errors. 



the analytical calculations of Section 2.2.2 and Fig. 5. In the 
calculations that have led to Figs. 3 and 4, the local beam of 
non-coherent curvature radiation is emitted, with the real- 
istic size appropriate for a preselected frequency z^, and for 
local p (as measured in the inertial observer's frame lOF). 
However, instead of using the bandwidth-integrated beam, 
in Figs. 3 and 4 we have used the 'bolometric beam' which 
is integrated over the entire frequency spectrum. This is be- 
cause: 1) the bolometric beam is given by the simple analyti- 
cal formula of eq. 8 in DRDIO, and 2) the shape of the bolo- 
metric beam happens to be similar to the observed shape 
of the BEC in J1012+5307, despite the latter is bandwidth- 
limited only. In this way we have avoided some deeply-nested 
calculation loops that would have otherwise been required 
to perform bandwidth-limited integration of the frequency- 
resolved beam, which cannot be expressed in a simple ana- 
lytical way. Such a realistic beam shape (ie. one which cor- 
responds to the observed bandwidth) is only used in Section 
3 to fit the double notches of B1929-fl0. The Lorentz fac- 
tor 7 in the bolometric beam is set to a value that corre- 
sponds to the realistic opening angle of bandwidth-limited 
beam at the central frequency of p — 820 MHz, and which 
also corresponds to the local lOF curvature radius of elec- 
tron trajectory. [Note: In the low-frequency limit, ie. when 
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the EEC's spectrum extends to much higher frequency than 
the observation frequency v, the frequency-resolved beam of 
curvature radiation does not depend on the Lorentz factor, 
so the observation frequency and the curvature of trajectory 
fully define the opening angle.] Thus, whereas the shape of 
the beam used for calculations illustrated in Figs. 3 and 4 is 
approximate, its opening angle is a realistic one. The aber- 
ration and propagation time effects (for fiat space-time) are 
included in these calculations. 

Note that on the LS of polar cap, the lOF radius of 
curvature increases with the radial distance r at a slower 
rate than the standard dipolar radius of curvature pdip- The 
latter is equal to: 



Pdip : 



(1) 



where a = sinOm/ sinOio is the fieldline footprint parameter, 
9m is the magnetic colatitude (measured from the dipole 
axis), 6io is the magnetic colatitude of the last open field 
lines at radial distance r, and i?ic is the light cylinder radius. 
The corresponding radius of curvature of electron trajectory 
in the lOF is: 

-1 

S _ . \ „ / -Rlc 




+ 2 sin Q = R] 



pdi] 



+ 2 sin a 



-1 



(2) 



This equation only refers to the leading side of polar tube 
and holds when the viewing angle C, is equal to the dipole 
inclination (a ~ C,). Comparison of eqs. ([ij and ([2} tells 
us that the increase of radius of curvature with altitude is 
smaller for the electron trajectory, than it is for the cor- 
responding B-field line along which the electron moves in 
the corotating frame CF. For the orthogonally-rotating MSP 
with P = 5.25 s (i?ic = 25 ■ 10^ cm), the dipolar radius of 
curvature at the leading side of polar cap is pdip = 6.7 ■ 10® 
cm, whereas the radius of curvature of the corresponding 
lOF trajectory is equal to p = 4.3 • 10® cm. In our numerical 
code, the lOF values of p are used. Within a limited part of 
magnetosphere they are given by eq. ^ (Dyks, Wright & 
Demorest 2010; Thomas & Gangadhara 2005). 

A sample result, presented in Fig. [3] has been obtained 
for an almost orthogonal rotator, with dipole tilt of a = 85°, 
C = 85°, P = 5.25 ms, v = 0.5, 0.8, and 1.4 GHz. The stream 
was infinitesimally thin and constrained to a single B-field 
line with the footprint parameter s = 1.2 and the surface 
magnetic azimuth 4'm = 95° . For the model of Fig. [2] to 
work, the radial profile of emissivity may be represented 
by any non-monotonic function of altitude, with a single 
maximum at some height. A convenient function that we had 
arbitrarily selected was the Gauss function of the lOF length 
of electron path (measured along their lOF trajectory). The 
emissivity had a maximum at some 'path-tracing' distance of 
^max from the polar cap, and the Icr width of the emissivity 
profile was equal to AZ = 0.1-Ric. A weak, residual RFM 
of r oc v~^'^^ was assumed, ie. the peak emissivity for the 
frequencies v — 0.5, 0.8 and 1.4 GHz occurred at distances 
of 0.31, 0.29, and 0.27i?ic from the polar cap. 

Results of the simulations shown in Fig. [3] confirm the 
expectations based on Fig. (2] because of the RFM, the ratio 
of fluxes in the two peaks becomes inverted, but the BEG 
stays fixed at the same pulse phase at all frequencies. The 
separation of maxima in the BEG increases for decreasing 
frequencies (Abfc oc i/"^^''), which is caused purely by intrin- 
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Figure 3. Frequency evolution of a BEC, calculated for a = = 
85°, </>m = 95°, s = 1.2, P = 5.26 ms. The extent of emissiv- 
ity along the electron path was Al = 0.1 -Ri^ and the RFM was 
r oc u^*^'^^. Note the asymmetry inversion, and the increase of A 
for decreasing w. This BEC is created by the viewing geometry 
represented by the lowest horizontal line in the next figure. 



sic properties of the emitted microbeam; it is not caused by 
the RFM. The trailing half of the simulated BEC is wider 
than the leading one (if both are measured at half max- 
ima). This results from the fact that the leading part of the 
BEC on average originates from a larger r than the trailing 
part (see Fig. [2]). In the simulation shown, the asymmetry 
of width amounts to ~ 4%, which is noticeably larger than 
observed for J1012. A close inspection of fig. 8 in DRDIO 
reveals asymmetry of width of ~ 1% magnitude. This may 
suggest that the BEC of J1012 originates from non-dipolar 
(more circular) B-field lines, or that the cut angle 5cut is 
larger than in the simulation (or both). 

The degree of asymmetry of peak flux a.t v Vsym, de- 
pends in a degenerate way on several factors, such as Al, 
Scut, and the strength of RFM. To measure the flux asym- 
metry, let us introduce the ratio Rf = Fi/Ft, where Fi and 
Ft denote the peak flux at the leading and trailing maxi- 
mum in the BEG. At frequencies far from fgym the ratio Rf 
can evolve with v or stay constant, depending on the radial 
profile of emissivity. This is because for a fixed viewing ge- 
ometry Rf simply represents the ratio of emissivities at two 
fixed altitudes: Rf{v) = F^{rA) /Fi,{rB), where the indices 
'A' and 'B' refer to the points in Fig. O 

For small values of Al or Scut the asymmetry Rf 
strongly depends on the viewing angle which can vary 
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as a result of precession, as well as on a. A change of by 
just 1° (from the value of 85° shown in Fig. O to ( = 84°) 
changes the flux ratio from Rf ~ 1 (Fig. [Sb) to Rf ~ 0.4. 
With the additional asymmetry in the width of the peaks, 
such components would probably be not identified as 'bifur- 
cated' at all (they look just 'double') . 

2.2 Apparent magnification by geometric effects 

The double features can be interpreted as a direct conse- 
quence of structured shape of microscopic emission beam. 
Let us consider that part of the curvature radiation beam 
which is polarised orthogonally to the plane of electron tra- 
jectory. Let the frequency of observation be much smaller 
than the peak frequency fcrv of standard curvature spec- 
trum: Uciv ~ 7 GHz 7^/p[cm], where 7 is the electron 
Lorentz factor. Radiation of this type has the intrinsic bi- 
furcation angle of: 

2^ = 0.8°{p7iy9)-'^\ (3) 

where p7 = p/(10^ cm), and fg = GHz). For dipolar 

curvature radii p7 ~ 1, the beam is up to 10 times nar- 
rower than the size of features observed in B0950, B1929, 
and J1012. If the physical model is correct, it implies that 
either p ~ 10* - 10^ cm, or that the observed width has been 
strongly enlarged by geometrical magnification. 




Figure 4. Sky-projected emission from 72 streams with s = 1.2 
and = 0°, 5°, 10°, 15°, etc. For all streams Zmax = 0.25iJic 
and Al = 0.1-Ric. The other parameters are the same as in the 
previous figure. Note the bifurcated nature of the emission stripes, 
and their tendency to assume horizontal direction at large \<j>\, 
that results in a small cut angle Scut- 



2.2.1 Probability of small- angle cut through the stream 

To illustrate the effects of geometrical enlargement, we have 
calculated a sky pattern of bifurcated emission from 72 
streams that emerge from the surface of neutron star at 
equal intervals of magnetic azimuth /S.<j)m ~ 5°. The re- 
sult is shown in Fig. [4l where the bifurcated nature of the 
emission can be discerned as a low-flux fissure at the cen- 
tre of each emission stripe. We have selected B-field lines 
with the footprint s = 1.2 to maximize the bifurcation an- 
gle (so it is more easily discernible in the figure) as well 
as to avoid some overlapping and pile-up of emission due 
to the caustic effects (which tend to occur for s<l). The 
most vertical emission stripes correspond to (prn ~ 190° (the 
top one), and (/f>m — 350° (bottom). The radial emissivity 
profile was the Gauss function with imax ~ 0.25i?ic and 
Al — 0.17?ic. The three horizontals mark sightline traverse 
at C = 20, 50, and 85°. The dipole tih was a = 85° so 
the radiation from the magnetic pole is located in Fig. [4] at 
(j} = ~2Rtis/Ric = -4.5°, C = 85°. 

As can be seen in Fig. 2] with increasing distance from 
the polar cap the emission of streams tends to assume hori- 
zontal direction and the stripes of emission are cut at small 
angles by the line of sight. Far away from the pole (|<^| ~ 90°) 
this is true for any viewing angle For ^ ~ a, the cut an- 
gle is small within the entire phase range shown in Fig. |4l 
except from the close vicinity of (j) — 0. For small ^ ~ 20° 
the observed bifurcation is large for any phase, because of 
the classical magnification by the 'not a great circle' effect. 
Thus, for the case shown in Fig. 31 it is only for moderate 
viewing angles ~ 30 — 60° , and for a limited range of phase 
101 < 70°, where the BECs are not strongly magnified. More- 
over, if we assume that easily-detectable pulsars are viewed 
at a small impact angle ((^ ~ q), then the strong enlarge- 
ment should be expected for a wide phase interval. This is 



because in this case the cut angle Scut is noticeably large 
only for limited phase range around (/> ~ 0, and (j) ~ 180° 
(see the dashed line for = 75° in Fig. [5] further below). 

We conclude that there is a fairly large parameter space 
in which the observed width of BECs can be considerably 
increased by geometric effects of either 1) the small cut angle 
Scut, or 2) the small viewing angle The presence of the 
first factor considerably increases the size of the space, in 
comparison to the case of conal beam, for which only the 
second type of enlargement was possible. 

2.2.2 Peak separation as a function of viewing geometry 

If the aberration and retardation are neglected, the appar- 
ent separation A of maxima in a BEC can be derived by 
application of spherical trigonometry to the standard dipo- 
lar magnetic field. The result is: 

cos K — cos^ ^ 



sin^ ( 



where 



sini/) 



(4) 



(5) 



sin Scut 

The bifurcation angle 2tf} is given by eq. ((3|, whereas the cut 
angle Scut is equal to: 

<5cut = 90° - Vpa, (6) 

where tppA is the usual polarisation angle between the sky- 
projected B-field and the sky-projected rotation axis: 



tantp-p 



(7) 



cos(<^ — (po) COS ^ sin a — cos a sin C 

(Radhakrishnan & Cooke 1969). Equation (|4]) is the split- 
fan analogue of the well-known equation for opening angle 
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Figure 5. Observed separation of peaks in a double feature pro- 
duced by tile sigiitline's traverse tlirougli a split-fan emission 
beam witii bifurcation of 2-)/> = 1°. Curves of different types are 
calculated from eq. Q with q = 80°, for six values of C = 20° 
(solid), 50° (dotted), 75° (dashed), 80° (dot-dashed), 85° (three- 
dot-dashed), and 90° (long-dashed). Note the obvious possibility 
of observing A > 5° within \4>\ < 120°, despite V = 0.5°. 

of conal beams. When the bifurcation angle 2tp, and the 
observed separation A are small, eq. Q reduces to: 



sin ( sin 5cut 

as was guessed in DRDIO. 

With the special relativistic (and altitude-sensitive) ef- 
fects ignored, the cut angle Scut is uniquely determined for 
every pulse phase as soon as a, ^ and (jio are selected. 
Therefore, eq. ^ actually gives the separation A as a func- 
tion of pulse phase 0. Fig.[5]presents examples of such curves 
calculated for a = 80°, tp = 0.5°, and various viewing an- 
gles (. It can be seen that for < a our sightline's traverse 
becomes instantaneously tangent to the projected B-field 
lines, which causes A to increase infinitely. For small ( <^ a 
this happens near 101 ~ 90° whereas for ( approaching a 
this phase interval moves closer to the dipole axis phase at 
= 0. 

In spite of the rather small intrinsic bifurcation of 2-0 = 
1°, the observed bifurcation A can reach large values within 
wide phase intervals that precede and lag the phase of dipole 
axis. If measured from the brightest (main) feature in a pulse 
profile, the known double features are located at phase lon- 
gitudes: 103.5° (B1929-hl0), 51°, 70° (J0437-4715), -30° 
(B0950+08), -55° and 130° (J1012-f 5307). As can be seen 
in Fig. O all these numbers, except from the last one, lay 
within the phase interval of strong magnification. Unfortu- 
nately, it is the last location which corresponds to the bright 
BEC of J1012-I-5307, the widest-observed so far. However, 
even this case can be accounted for by the model, if the 
dipole inclination a in J1012-I-5307 is very small or close to 
90° (in the latter case the phase (p — 130° becomes equiv- 
alent to <j) — —50°, and there is no difference between the 
main pulse and interpulse). 

2.2.3 Possibility of interpulses 

Regions of strong radio emission in pulsar magnetosphere 
are definitely more localised than shown in Fig. [l] because 



the observed profiles usually have duty cycles much smaller 
than 360°. However, most of the emission that contains dou- 
ble features is not localized within a narrow range of phase. 
This is the case of B1929+10, B0950+08, and J0437-4715. 
Narrow radio pulses become a possibility when azimuthal 
limitations for the extent of emission are allowed, or when 
the radial extent of strong emissivity is not too large, yet 
noticeable, say Ar <0.1-Ric. 

It is worth to note that even for very extended emission 
regions, interpulses do not have to be always detectable, 
because of bulk topological properties of pulsar magneto- 
sphere. One example of this is the outer gap model, in which 
'the other pole' emission is undetectable by a single observer 
because it is constrained to regions above the null-charge 
surface (Romani & Yadigaroglu 1995; Wang et al. 2011). An- 
other example is the slot gap model (Arons 1983; Muslimov 
& Harding 2003) with the large-scale acceleration possible 
only on favourable (poleward) magnetic field lines. Thus, in 
a large fraction of existing magnetosphere models the invisi- 
bility of the other pole is an inherent property. It is therefore 
possible that radially-extended emission streams are present 
in many objects, but are not directly manifested by unavoid- 
able second-pole emission. 



3 PIECE- WISE FITS OF PHYSICAL CURVES 
TO THE DOUBLE NOTCHES 

As can be seen in Fig. [21 the altitude of detectable emission 
changes gradually while the sightline is obliquely traversing 
through the double pattern. Therefore, the width and height 
of the sampled microbeam also gradually change with pulse 
phase, because they depend on the curvature radius p(r) and 
r-dependent emissivity. Obviously, the same is true for dou- 
ble features in the absorption version, ie. for double notches 
which may have r-dependent width and depth. Precise fits 
to double features then require a full 3D modelling of the 
stream-cut geometry, ie. the r-dependence of the depth of 
the notches, and the r-dependence of p need to be assumed 
in addition to q, and the specific functional form of the mi- 
crobeam. This kind of calculation is feasible with the code of 
type that produced the result shown in Fig. [4] However, the 
frequency-integrated (bolometric) microbeam which is now 
used by the code, would need to be replaced by an observing- 
band-integrated beam. And the code would have to be em- 
bedded into a fitting routine with some iterative capability 
of heading towards minimum. Instead of these numeri- 
cal developments, we make a direct, but piece-wise fit to the 
double notches of B1929+10 using the bandwidth-integrated 
beam. More precisely, it is assumed that the left half of the 
notches on average originates from some unspecified altitude 
with one curvature radius p{rA), and the other half on av- 
erage originates from another altitude and pirn). Since the 
scale of the fitted beam depends only on the curvature ra- 
dius p, it is only the values of pi that are provided by the 
fit. Any further association of pi with radial distances r^ is 
not unique because of unknown viewing geometry (5cut). As 
will be seen below, the limited quality of the best-available 
data on the notches of B1929-I-10, makes this procedure rea- 
sonable, because fractional values of /dof can be reached 
even with this simplified model. 
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3.1 Functional form of the fitted microbeam 

The radio flux at pulse longitudes surrounding the notches of 
B1929+10 changes linearly (see fig. 1 in DRR07). Therefore, 
before making the fit we remove this linear trend to get a 
constant flux normalised to unity (Fig. [51). Then we use the 
following function to fit the data: 



(9) 



where the integration is over the observed frequency band- 
width and 77crv represents the shape of curvature radiation 
microbeam: 



^11 + Vl 



+ 



where 



^ = 1/7^ + 7/)^, and 
3c 



(10) 
(11) 
(12) 

(13) 
(14) 



(eg. Jackson 1975) . The symbol p denotes the radius of cur- 
vature of electron trajectory in the lOF, uj — 2-kv, c is the 
speed of light, and K's are the modified Bessel functions. 
The angle is the azimuth of the line of sight measured 
in a frame with z-axis along the instantaneous acceleration 
vector (it is the azimuth of sightline measured around the 
vector of curvature radius) The two terms in the square 
bracket represent decomposition of the beam into two parts: 
one polarised parallel, and the other polarised orthogonally 
to the plane of electron trajectory. In the absence of the par- 
allel mode the normalisation constant Cy is equal to zero, 
and the microbeam has the double-lobed shape with zero 
flux in the centre. Therefore, if spatial convolution effects 
are not included in the flt, the central maximum of the dou- 
ble notches, which does not reach the 'off-notch' level, is 
not well reproduced (Fig. [5^, with x^/dof = 1.1 for the left 
notch, and /dof = 0.6 for the right notch). The decreased 
flux at the central maximum probably results from some 
convolution effects, such as the non-zero width of the non- 
emitting/eclipsing region. To avoid the complexity of such 
convolution (eg. the selection of transverse plasma density 
profile), we instead lower down the maximum by adding a 
bit of the parallel mode (Cy ~ 0.13Cx, Fig.|6jD), which has a 
single-peaked pencil-like shape. Note that the data may well 
contain no parallel mode at all. We simply use the non-zero 
value of C|| as a convenient way to improve the fit. The spa- 
tial convolution seems to be the most natural explanation 
for the low flux at the center of 'W. In the fltting proce- 
dure we simply set ip = (j> ~ <j)c, where (f)c is the phase of the 

^ Let us introduce a plane that is orthogonal to the plane of 
electron trajectory, and that includes the instantaneous velocity 
vector V. Let us project the line of sight onto this orthogonal 
plane. ij) is the angle betweeen the projection and the plane of 
electron trajectory. Note that -0 is NOT the angle between the 
trajectory plane and the sightline itself, as was wrongly declared 
in DRDIO. 



V 



V 



-20 -16 -12 




-20 -16 -12 



Figure 6. A fit of classical curvature radiation beam (eq. [9J to 
double notches in the 327-MHz pulse profile of B1929-I-10. Each 
notch is fitted separately, within the left-hand side, and right- 
hand side pair of solid vertical lines. The pair of dotted verticals 
marks the phase interval used to remove the linear increase of 
flux with phase. In panel a) the orthogonal polarisation mode is 
only used (Cy = 0). In panel b) we have fixed Cy ~ 0.13 to 
match the central data points. The data were collected with the 
Arecibo Telescope (DRR07); the bandwidth was 25 MHz. 



centre of the 'W. In the low- frequency limit oi v <^ Vav, 
the shape and size of curvature radiation beam do not de- 
pend on the Lorentz factor 7 of emitting particles. Since 
the frequency of observation shown in Fig. [6] is relatively 
low {f — 327 MHz), it is reasonable to perform the fitting 
without any reference to 7. To minimize the number of fit 
parameters in this way, we simply set 7 to some 'very high' 
value (eg. 10^) to ensure that v ^ i'crv(7,p). The width 
scale of the beam is then changed only through changes of 
p. Therefore, the fitted value of curvature radius also incor- 
porates any influence of viewing geometry on the scale of the 
notches. Thus, the fit actually provides us with the value of 
Pmin ~ p sin'' ^ sin'^ (5cut , where the subscript 'min' has the 
meaning of minimum value of p for both ^ and 5cut equal 
to 90"^. In spite of the wide bandwidth {Au/v = 8%), the 
frequency integration in eq. ([9)l does not change the beam 
shape considerably. The bandwidth-integrated beam stays 
similar to the frequency-resolved shape of eq. (|12p . 

As can be seen in Fig. |6l the curvature radiation beam 
reproduces the notches of B1929-I-10 very well. Even when 
no account is done for the lowered flux of the central mini- 
mum (Fig. [6^), the shape of the curvature subbeam neatly 
reflects the steepness of the outer walls of the 'W, as well 
as the bulk width of the minima, as compared to the sep- 
aration between them. With the central maximum lowered 
down for the price of one additional fit parameter (Cy), the 
quality of the flt reaches x^ /dof = 0.7 and 0.2 for the left 
and right-hand side notch respectively (Fig. (BJs). 
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3.1.1 Curvature radii 

The fitted values of curvature radius are: pmin = 1.17, and 
0.85 • 10^ cm for tlie leading and trailing notch, respectively. 
These values assume orthogonal and equatorial cut through 
the stream (5cut = C = 90° ) and should be understood as the 
lower limits for p = pmin/(sin^ C,s\v? Scut). The value of pmin 
fitted for the leading notch is somewhat larger than for the 
trailing one. This seems to be natural for the double notches 
of B1929+10, because small (5cut is likely for this object: the 
interpulse suggests that a ^ C, ^ 90°, and the notches are 
located 103.5° after the main pulse, ie. within the region of 
strong amplification through the small-5cut effect (Fig. [S]). 
Note that already for moderately small 5cut = 28° (and 
1^ = 90°), the inferred value of p reaches lO'^ cm, whereas 
for (5cut = 13°, p — lO' cm. The fitted size of the double 
notches (A = 5.37° ai v = 327 MHz) is then compatible 
with dipolar values of near-surface p in the closed field line 
region. 

In the plane of the magnetic equator the curvature ra- 
dius of dipolar B-field lines is equal to: 

p,, = p{e^ = 90°) = ^, (15) 

where 6m is the magnetic colatitude of the emission point, 
measured from the dipole axis. Thus, for every pulsar with 
radius R^s = 10^ cm, the minimum curvature radius avail- 
able in its dipolar magnetosphere is equal to p = 3.3 • 10^ 
cm. 

In B1929-I-10 the double notches are observed at ~ 
103.5°, ie. roughly at the right angle with respect to the 
dipole axis. Therefore, it is useful to find the surface value 
of p at the position where the B-field is pointing at the right 
angle 9b ~ 90° with respect to the dipole axis. The mag- 
netic colatitude for this point is equal to 8m — 54.73° (more 
precisely: = 2/3), and the corresponding curvature 

radius is: 

px = p(eB = 90°) = ^r, (16) 

which at the star surface gives p — 0.866-Rns — 10^ cm. 
Thus, even in the presence of only modest geometrical am- 
plification ((Scut — 30°, ^ ~ 90°) the dipolar curvature radii 
in the closed field line region are consistent with the observed 
width of double notches in B1929+10. 

In the unlikely case of large Scut, the substellar values of 
pmin would have to be interpreted in terms of small-scale dis- 
tortions of dipolar B-field (Harding & Muslimov 2011; Ru- 
derman 1991). Another option, at least in principle, would 
be the gyration radius of primary particles near the light 
cylinder: 

77? ^ 

rs = 1702 cm (17) 

nie B 

where m/me = m' is the mass of radiating charges in units 
of the electron mass m^, and B is in Gauss. This radius 
translates into the bifurcation angle of: 

2V = 14.4°f4^y'', (18) 

which is broadly consistent with observations for B < 10^ 
G and 7> 10®. However, to produce strong and unsmeared 
double features, the gyrational motion would have to be spa- 
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Figure 7. Average merging-rate index a in the relation A oc u~'^. 
'BEC refers to bifurcated emission components, whereas 'DN' to 
double notches. Dashed verticals mark the range of a expected if 
the doubleness is caused by the double-lobed microscopic beam 
of curvature radiation (eq. The notches of B0950-I-09 (not 
included in this figure) are also consistent with a = 1/2 (see fig. 9 
in DRR07). 

tially ordered within a large volume of magnetosphere. This 
does not immediately seem to be natural. 

3.2 Separation as a function of frequency 

In view of the remarkable agreement seen in Fig. [B] it is 
needed to verify if the separation A between the minima (or 
maxima) in double features is also consistent with the prop- 
erties of curvature radiation. According to the classical elec- 
trodynamics, at low frequencies {f <^ fcrv) the separation 
A decreases with increasing u according to eq. Q, ie. ac- 
cording to the power-law A oc i^"" with the index a — 1/3. 
The bifurcation angle in this case is i/i 3> I/7, where 7 is the 
Lorentz factor of the emitting particles. At the high-energy 
end of the emitted spectrum, ie. within the high-energy 
spectral decline, the features merge faster, with the index 
a = 1/2 (eg. Jackson 1975), and the bifurcation becomes 
very small: ip ^ I/7. A simple arithmetic average then al- 
lows us to expect an 'average' merging rate of a = 0.417 
around the spectral maximum. It is worth to emphasize that 
many models of coherent curvature emission follow the same 
tp{v) as the non-coherent curvature radiation. Such models 
include those based on bunching of charges (Ruderman & 
Sutherland 1975), as well as some maser mechanisms (Luo 
& Melrose 1992; Kaganovich & Lyubarsky 2010), 

The theoretical range of merging rate for the curvature 
radiation is compared to the available data on double fea- 
tures in Fig. [T] which includes data on both the emission 
(BEC) and absorption features (DN). The figure takes into 
account the new data on J0437— 4715 from Yan et al. (2011), 
as well as unpublished low- frequency data on this object. 
For this reason the two points for J0437— 4715 lay closer to 
a = 1/2 than in fig. 6 in DRDIO. 

One can see that within the statistical Icr errors, the 
measured merging rate is consistent with the one expected 
for the curvature beam. Inspection of A as a function of u 
may even be revealing some decrease of a at lower frequen- 
cies (see fig. 6a in DRR07). We remind here that it is not the 
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merging rate that differentiates the curvature radiation from 
the other emission processes, because the latter can also have 
a<l/2 (eg. Melrose 1978; DRDIO). However, ff the origin 
of the observed doubleness is not related to the microscopic 
beam, this would possibly be manifested by merging index 
values from outside of the range [1/3, 1/2]. The large num- 
ber of indices near a = 1/2 may be considered to be not 
quite natural for the curvature radiation, because double 
features should be more easily observable at low i/, where 
the beam is wider, and the merging index a is closer to 1/3. 
However, double features at high v can be more easily af- 
fected by blurring effects, which can also bring the maxima 
increasingly closer together, thus increasing a. We therefore 
conclude that the analysis of merging rate is broadly consis- 
tent with the curvature origin of double features. 

3.3 Fitting the BEC of J1012+530T 

We have also made several fits to the BEC of J1012-f5307, 
in which the outer wings are much less steep than the outer 
wings of the double notches in B1929-I-10. Therefore, we 
have used functions that are a convolution of some bell- 
shaped plasma-density profile (like the Gauss or Lorentz 
functions) with the microbeam of eq. ([9]). This type of model 
has appeared to be unable to reproduce the BEC: before the 
flux in the wide outer wings has been reached, the flux in the 
central minimum was already too large. We conclude that 
with the use of the microbeam from eq. ((9]l, the shape of 
the outer wings of the BEC in J1012-I-5307 cannot be repro- 
duced by convolving it with bell-shaped density functions. 
Another convolution that is possible for the simple case of 
f <^ i^cTv, and which remains to be tried out, is a spread 
in curvature radii p, which can act similarly to the inte- 
gration over a wide frequency range (much wider than the 
bandwidth). Yet another possibility is the case of v ~ Vav, 
ie. the observation frequency being close to the characteristic 
frequency of curvature spectrum, in which case the opening 
angle is V' ~ 1/7, where 7 is the electron energy in units of 
mc^. Finding the shape of BEC in such a case would require 
us to make integration over the electron energy distribution 
(in addition to the spatial convolution). 



4 CONCLUSIONS 

We conclude that the stream-cut geometry is a successful 
model for such components in radio pulse profiles like the 
apparently 'conal' components in J0437— 4715, or the BEC 
in J1012-I-5307. Such geometry naturally explains the fixed, 
frequency-independent location of these features in the pro- 
files. The components stay at a fixed phase not because of 
the 'lack of REM', but because the emission region is nar- 
row in azimuth (stream-like emitting region, or spoke-like, if 
there are more streams). If the REM is added to this picture, 
it does not move the components away from the profile cen- 
ter. Instead, in the case of non-orthogonal cut through the 
stream, the REM is making the double features asymmetric 
at all frequencies, except from at i^sym. This conjecture is 
valid for stream-like emitters with REM regardless of the 
actual origin of the doubleness, ie. regardless of whether the 
bifurcation has the micro- or macroscopic origin. 

If we assume that the doubleness refiects the shape of 



the microphysical beam, then the curvature radiation in the 
extraordinary mode becomes a natural process, because: 1. 
it is intrinsically double-lobed (bifurcated), which is a neces- 
sary condition for producing deep double notches (DRDIO); 
2. the extraordinary part of the curvature beam reproduces 
the shape of notches in B1929-I-10; 3. the rate of merging of 
double features with increasing frequency is consistent with 
the curvature origin; 4. the width-scale of double features in 
B1929-H0 and J0437-4715 is consistent with dipolar cur- 
vature radii; 5. the high polarisation of double features is 
consistent with the single-mode nature of the double-lobed 
beam. 

This said, it is not quite clear why some of the observed 
double features have the width of > 7° at 1 GHz, despite the 
intrinsic beam is expected to be only 1° wide. However, the 
quite large probability of strong geometric magnification, 
along with the possibility of smaller-than-dipolar p, do not 
really allow us to consider this issue as a serious problem for 
the curvature microbeam model. As recently noted by Yan 
et al. (2011), who verify fairly large sample of high-quality 
profiles of MSPs, bifurcated features are rare among MSPs. 
This is in line with the special requirements that need to 
be fulfilled for these features to be observable. The require- 
ments apparently include at least the following: 1. the spa- 
tial extent (azimuthal width of the stream) must not be too 
large to blur the feature; 2. the feature be best magnified by 
geometrical effects to be wide and easily detectable. It is not 
clear why we do not observe features with A< 1°, that ful- 
fill the condition no. 1, but are not magnified geometrically. 
Generally, however, the intrinsic narrowness of the curva- 
ture beam seems to be consistent with the rarity of double 
features in the observed profiles. This is because for even 
a small azimuthal extent of emission region, the doubleness 
should become hidden by spatial convolution effects. 

We conclude that the asymmetry of double features at 
different frequencies, the shape of notches in B1920-I-10, the 
way they merge, their high polarisation, and their fixed 
phase locations can all be understood in terms of sight- 
line cuts through streams of plasma that emits curvature- 
radiation in the orthogonal polarisation mode. 
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